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ABSTRACT 

We present data from simultaneous Chandra, XMM-Newton and BeppoSAX 
observations of the Seyfert 1 galaxy NGC 3516, taken during 2001 April and 
November. We have investigated the nature of the very flat observed X-ray 
spectrum. Chandra grating data show the presence of X-ray absorption lines, 
revealing two distinct components of the absorbing gas, one which is consistent 
with our previous model of the UV/X-ray absorber while the other, which is 
outflowing at a velocity of ~ 1100 km s -1 , has a larger column density and is 
much more highly ionized. The broad-band spectral characteristics of the X-ray 
continuum observed with XMM during 2001 April, reveal the presence of a third 
layer of absorption consisting of a very large column (« 2.5 x 10 23 cm -2 ) of highly 
ionized gas with a covering fraction ~ 50%. This low covering fraction suggests 
that the absorber lies within a few lt-days of the X-ray source and/or is filamen- 
tary in structure. Interestingly, these absorbers are not in thermal equilibrium 
with one another. The two new components are too highly ionized to be radia- 
tively accelerated, which we suggest is evidence for a hydromagnetic origin for the 
outflow. Applying our model to the November dataset, we can account for the 
spectral variability primarily by a drop in the ionization states of the absorbers, 
as expected by the change in the continuum flux. When this complex absorp- 
tion is accounted for we find the underlying continuum to be typical of Seyfert 
1 galaxies. The spectral curvature attributed to the high column absorber, in 
turn, reduces estimates of the flux and extent of any broad Fe emission line from 
the accretion disk. 
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Subject headings: galaxies: active - galaxies: individual (NGC 3516) - galaxies: 
nuclei - galaxies: Seyfert 


1. Introduction 

It is now understood that intrinsic UV and X-ray absorption is present in the spectra 
of at least half of all Seyfert 1 galaxies (Reynolds 1997; George et al. 1998; Crenshaw et 
al. 1999). High resolution UV spectra obtained with the spectrographs aboard the Hubble 
Space Telescope and, most recently, X-ray spectra obtained with Chandra and XMM have 
revealed blueshifted absorption lines, signifying massive outflows of gas from the active nuclei 
of these galaxies (Crenshaw, Kraemer, & George (2003), and references therein). X-ray and 
UV spectra of NGC 3516 (z=0. 008836; Keel 1996) have shown evidence for a significant 
column of ionized gas along the line-of-sight to the nucleus (Kriss et al. 1996; Netzer et al. 
2002; Kraemer et al. 2002). Recent observations have found NGC 3516 to be in a relatively 
low flux state. Chandra LETG data suggested a model whereby a constant column of gas 
with Nh ~ 8 x 10 21 atom cm -2 reacts to changes in the nuclear ionizing flux (Netzer et al. 
2002). Kraemer et al. (2002) argued that the strong UV absorption lines, e.g., Lya, N V AA 
1238.8, 1242.8, and C IV AA 1548.2, 1550.8, arose in the same column of gas. 

We present here Chandra HETG spectra, XMM-Newton (hereafter XMM) EPIC CCD 
spectra and PDS data from simultaneous BeppoSAX observations of NGC 3516 from 2001 
April and November. Previous analysis of a subset of these data concentrated on just the 
Fe Ka regime in data from 2001 November. The Fe K a line showed interesting structure 
and rapid evolution during the 2001 November observation, and the combined HETG/EPIC 
data in the hard X-ray regime have been detailed by Turner et al. (2002). Initial analysis of 
the 2001 November data also revealed a very flat observed continuum shape. November data 
from the HETG showed no significant absorption lines, leaving ambiguity as to the nature 
of the flat observed spectrum that could be either intrinsically flat; have a ’normal’ photon 
index and appear flat due to complex layers of absorption or be dominated by reprocessed 
components mimicing a flat spectrum (i.e. a very broad line from the accretion disk (hereafter 
denoted a ’diskline’) of large equivalent width and associated Compton hump). As the 
ambiguity could not be broken using 2001 November data alone, Turner et al. (2002) used 
a simple flat powerlaw parameterization of those data, against which to examine the newly- 
discovered narrow and redshifted line emission from Fe. 

Analysis of data from 2001 April has now shown footprints of highly-ionized absorbing 
gas, never seen before in this source. We present a detailed parameterization of NGC 3516 
based on data from both April and November epochs of 2001. 
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2. The Multi-satellite Observations 
2.1. The Chandra Data 

Chandra observations of NGC 3516 were performed at two epochs during 2001 using the 
High-Energy Transmission Grating Spectrometer (HETGS; e.g. Marshall, Dewey & Ishibashi 
2004) in conjunction with the ‘S-array’ of the Advanced CCD Imaging Spectrometer (ACIS; 
e.g. Garmire et al 2003). We note that to-date, these are the only Chandra observations 
of NGC 3516 to have been performed with the HETGS/ACIS-S instrument configuration. 
The first observation (obsids 2431 and 2080) was performed over the period April 09 14:12 
to April 11 14:53 UTC, but with an interruption (April 10 01:17-17:55 UTC) when solar 
activity led to sufficiently high background levels that the detectors were turned off. The 
second observation (obsid 2482) was performed between 2001 November 11 01:00 UTC to 
November 12 02:19 UTC. From the times of the observation listed above, it can be seen 
that the second period of Chandra observations in April (~21 hr during obsid 2080) is 
simultaneous with XMM, as is ~10 hr of the observation in November. 

The Chandra data from both epochs were processed using the ciao (v3.0.1) and ftools 
(v5.3.1) software packages, along with calibration files from the CALDB (v2.26) release 7 . We 
also made use of our own software for a number of the analysis tasks. The data were 
processed in the standard manner, including the removal of bad detector pixels and other 
detector artifacts, removal of events with detector ‘grades’ not equal to 0, 2, 3, 4, or 6. 
The selection criteria used resulted in exposures of ~36 and 73 ks for the first and second 
segments of the April data, and 88 ks for the November observation. 

The background was negligible and so no background subtraction was performed for 
the HETGS spectra. The mean count rates in the 0. 7-8.0 keV band obtained for summed 
(positive and negative) l st -orders of the Medium-Energy Grating (MEG) arm were 0.110 ± 
0.023 and 0.211±0.027 ct s -1 for the two segments of the April observations, and 0.123±0.013 
ct s _1 during the November observation. The corresponding values for the High-Energy 
Grating (HEG) arm were 0.069 ±0.013, 0. 127 ±0.018, and 0.076 ±0.010 ct s -1 (respectively). 
In the spectral analysis presented below we consider the data from obsid 2080, as this is 
simultaneous with the good subset of the PN observation. 


including the ac isD 1999-08- 13contamN0003. fits file to correct for the time-dependent contamination 
on the surface of the CCDs (e.g. see Marshall et al. 2004). 



2.2. The XMM Data 


XMM observations of NGC 3516 were performed covering 2001 April 10 11:14 - April 
11 23:02 UTC and November 09 23:13 - November 11 10:54 UTC. EPIC data utilized the 
thin filter with PN Prime Small Window mode. These data were processed using SAS 
v5.3.3. and were screened to remove hot and bad pixels and periods of high background. 
For the April epoch we took only the subset of PN data which were simultaneous with the 
second piece of Chandra coverage (April 10 17:56 - 14:53 UTC), yielding a 49 ks exposure. 
This is the period during which EPIC background levels were lowest, and data outside of 
this time range resulted in data with inadequate subtraction of the background spectrum. 
In November data screening yielded an exposure of 81 ks. At the flux level of NGC 3516 
photon pileup is negligible in EPIC data. During the analysis we found the MOS spectra to 
have much lower signal-to-noise ratio than the PN data, thus we performed our analysis on 
the PN spectra alone. For the PN, instrument patterns 0-4 were selected and spectra were 
extracted from a cell ~ 0.94' diameter, centered on the source. Background spectra were 
extracted from a nearby region for the PN (5 — 10% of the source count rate). 

NGC 3516 yielded 5.48±0.01 ct s _1 during 2001 April and ~ 3.25 ± 0.01 ct s -1 during 
2001 November, over the 0.35 - 10 keV band. Table 1 shows source fluxes for the two epochs, 
and for the earlier LETG observation which we pulled from the archives and compared to 
our best-fitting model for extraction of fluxes. NGC 3516 occupied the lowest known X-ray 
flux state during the 2000 October and 2001 November epochs, and during the early part of 
the 2001 April Chandra observations (obsid2431), rising to a higher flux later in 2001 April. 
This time variability allowed us to track the source behavior with changing flux. 


2.3. The BeppoSAX Data 

BeppoSAX observed NGC 3516 2001 April 10 UT 16:05 - April 12 22:48 UTC and 2001 
November 09 UT 20:27 - November 12 08:01 UTC. We utilized data from the Phoswich 
Detector System (PDS; Frontera et al. 1997) providing spectral data over ~ 15 - 200 keV. We 
used products from the BeppoSAX database, created using a standard data reduction using 
the SAXDAS software version 2.0. Data from the four PDS phoswich units were combined, 
after gain equalization and linearization. The instruments were switched off during SAA 
passages and PDS data acquired during the first 5 minutes after each SAA passage were 
rejected to exclude periods of highly variable risetime threshold. Standard screening of the 
BeppoSAX data yielded events files with effective exposure times of 40 ks (April) and 63 ks 
(November). 
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The PDS had a collimated fov with FWHM = 1.4 degrees. The region around NGC 3516 
is free of contaminating hard X-ray sources of significant flux within this radius. The PDS 
instrument was rocked back and forth to obtain source and background data covering the pe- 
riod of the observation. The PDS rocking mode generally provides a very reliable background 
subtraction. We obtained, after background subtraction, 0.86±0.04 PDS cts -1 during 2001 
April; 0.68±0.03 PDS cts _1 during 2001 November. 


3. X-ray Spectral Analysis 

3.1. Chandra Evidence for Multiple Components of Absorption 

The HETGS data were compared to an absorbed continuum model which parameterized 
the curvature of the source, and then examined for residual features. In Fig. 1 we show some 
of the most interesting sections of the spectrum obtained during 2001 April, where we have 
combined the (flux-calibrated) ±l st -order spectra from the MEG and HEG arms and plotted 
the data using O.OlA bins (corresponding to the ~ la width of the spectral resolution of the 
MEG, and twice that value for the HEG). Immediately evident are several deep absorption 
lines from highly ionized species of several elements, most notably the H-like species of Mg, 
Si and S. 

Based on an analysis of a Chandra Low-Energy Transmission Grating Spectrometer 
(LETGS) spectrum of NGC 3516 obtained on 2000 October 6, Netzer et al. (2002) argued for 
the presence of a relatively low-ionization (log U = —0.72) absorber, with a column density 
log N h = 21.9. Kraemer et al. (2002) suggested the same gas was the source of the strong 
UV absorption lines of H I, C iv, and N v detected in contemporaneous HST/STIS medium- 
resolution echelle spectra. Since the 2001 April and November observations exhibited similar 
soft X-ray spectral characteristics, it is plausible that a column of low ionization gas (which 
we will refer to as the “UV component”) was still present. Assuming a sufficiently large 
velocity dispersion within the absorbing gas (we assume 300 km s -1 here), the ionization 
levels in such a component are sufficient to give rise to detectable absorption lines (due 
to inner-shell transitions) from Mg vn, vm, ix and Si vn, vm, ix. However, this UV 
component is not sufficiently-ionized to contain any H- and He-like species of these elements. 
Nevertheless, as can be seen from Fig. 1, we do see clear evidence for absorption lines from 
H-like Mg, Si and S. The presence of absorption lines from the He-like species is less clear, 
although the data are highly suggestive of a line due to Si xm, and possibly also Mg xi. 
Thus we have strong evidence of at least one other, much more highly-ionized, component 
of absorption. 
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In Fig. 2 we show the kinematical constraints that can be obtained from the Mg xn, 
Si XIII, Si xiv, and S xvi lines (the He-like lines of Mg and S provided no useful constraints, 
primarily due to the signal-to-noise limitations). For this analysis we adopted the rest- 
wavelengths quoted in Verner, Verner & Ferland (1996),Verner et al. (1996), which are in 
good agreement with other compilations. Specifically for the H-like Is 2 S \/2 — > ^P 2 Py 2 , 1/2 
doublets (with oscillator strengths in the ratio 2:1) we assumed wavelengths of 4.7274 & 
4.7328A for S xvi, 6.1804 & 6.1858A for Si Xiv, and 8.4192 & 8.4246A for Mg XII. For 
the Si XIII ls 21 ^ — >• ls.2p 1 P° transition we adopted 6.6480A. For each ion, we constructed 
high-resolution ’opacity spectra’ for various combinations of ionic column density, N( ion), 
and velocity dispersion, v a . Both absorption-line (assuming full Voigt profiles) and bound- 
free absorption were included in these grids, and we assumed that the absorbing material 
completely covers the underlying X-ray source (i.e. a covering-factor of unity). The grids 
were then used to construct model spectra by assuming an underlying continuum (see below). 
These spectra were then convolved using the instrumental spectral response and compared 
to the data using xspec (v.11.3.1), with the redshift of the opacity spectra (and hence the 
offset velocity, v 0 ff, of the absorbing material with respect that expected given the recessional 
velocity of the host galaxy). 

For the analysis presented in Fig. 2, for each ion we restricted the analysis to the data 
within ±10 4 km s _1 of the expected location (after correcting for the systemic velocity) of the 
Is — » 2 p transition, and to within the same velocity range of the corresponding bound-free 
edge. Given the limited signal-to-noise of the data in the vicinity of the Is — > 2 p lines, the 
data in the vicinity of the bound-free edge had to be included to prevent solutions in the low 
V(j , very-high N{ ion) region of parameter-space. The MEG and HEG data within the above 
spectral ranges were fitted simultaneously (with the different instrumental responses taken 
into account). We adopted spectral bins of width 0.010 and 0.005 A for the MEG and HEG 
(respectively) and used a C-statistic minimization technique within xspec. In all cases a 
simple, smooth, underlying continuum was assumed, the normalization of which was allowed 
to vary during the spectral fitting. 

The dashed curves in Fig. 2 show the 90% confidence regions in v a -v 0 ff space obtained 
from the individual Mg XII, Si xm, Si xiv and S xvi ions. It can be seen that the results 
for the individual Mg xn, Si xm and Si xiv ions are all consistent with v 0 ff in the range 
—400 < v 0 ff < —1400 km s -1 , but v a cannot be constrained. The results for S xvi appear to 
be slightly discrepant with —1450 < v 0 ff < —2050 km s -1 , and v a < 400 km s _1 . This dis- 
crepancy may be partly the result of statistical fluctuations in low signal-to-noise data or may 
be indicative of additional spectral complexity caused by additional kinematic components 
(see below). From a combined analysis of all four ions we found v 0 ff = —11401 3^ km s" 1 . 
By comparison, the UV absorption components show velocities up to ~ 1300 km/s (Kraemer 
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et al. 2002). 

In contrast, the mean outflow velocity from the lower ionization S, Si and Mg lines was 
close to systemic, consistent with an origin within the UV absorbers. The strength of Si XI 
line (see Fig. 1) may indicate the presence of an additional, distinct layer of gas. However, 
the relative weakness of Six and Sixii features led us to question the significance of the Sixi 
absorption line and thus we note it but did not attempt to model it. 

While Chandra grating spectra from 2001 April allowed detection of the discreet ab- 
sorption lines, enabling us to pick out previously-unresolved layers of gas, there was still 
some ambiguity as to the physical conditions (i.e, column density and ionization) of each 
layer. The lines were not resolved by HETGS and thus we needed to fit the line profiles to 
obtain further constraints on the gas. There is a significant trade-off between the assumed 
turbulent velocity and the column density required to produce each line, as discussed in the 
next section. 

Unfortunately the source was fainter during November and HETGS data from the 2001 
November epoch showed no significant narrow absorption features although those data are 
consistent with the presence of the absorption lines found in 2001 April. 


4. Modeling the Absorbers 
4.1. Inputs to the Models 

Our general approach to modeling is to start with a single-component absorption-model, 
including additional layers only as necessary. Hence, although it is entirely likely that the 
intrinsic absorption is more complex, we initially restricted our models to include only the 
two components whose footprints are clearly seen in the HETG spectra. Our next step was 
to generate model tables of the two zones of X-ray absorbing gas, and then we used these to 
fit the moderate resolution EPIC PN spectra, using the good signal-to-noise available there 
to constrain the broad spectral components. 

The photoionization models for this study were generated using the Beta 5 version of 
Cloudy (G. Ferland 2003, private communication). We modeled the absorbers as single-zoned 
slabs of atomic gas, irradiated by the central source. We assumed a similar spectral energy 
distribution of the continuum radiation emitted by the central engine to that described in 
Kraemer et al. (2002), with the exception that we modified the spectral energy index, a, 
above 500 eV based on our initial fitting of the 2001 April data, such that: a = —1.0 for 
< 100 eV, a = —2.3 over the range 100 eV < h.u < 500 eV, and a = —1.0 above 500 eV. 
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The models are parameterized in terms of the dimensionless ionization parameter U , which 
is the ratio of ionizing photons per H atom at the illuminated face of the slab. Based on 
our models of the X-ray emission lines in NGC 3516 (Turner et al. 2003), we assumed solar 
elemental abundances (e.g. Grevesse & Anders 1989), with the exception of nitrogen, for 
which we adopted an N/H ratio of 2.5 times the solar value. The absorbing gas is assumed 
to be free of cosmic dust. 

We initially generated a grid of photoionization models around values of U and N #- 
suggested by the Chandra detections of the H- and He-like lines discussed above (which we 
will refer to as the “High Ionization” component). For simplicity, we assumed a 300 km s' 1 
velocity dispersion for the model, consistent with HETG constraints on the gas, and applied 
the predicted ionic column densities to fit the absorption-line profiles, using the wavelengths 
and oscillator strengths from Behar & Netzer (2002) and Verner, Verner & Ferland (1996). 
The predicted Sixiv and MgXII lines from our “best-fit” model (logl/ = 1.7, Log Ng = 22.2) 
were compared to the Chandra spectra and showed reasonable agreement between the model 
and the data. Velocity spectra (Fig. 3) indicate a mean bulk velocity -1100 km s' 1 for the 
High Ionization gas component. 

While the lower-ionization inner shell lines were expected to arise in the UV absorber, 
the signal-to-noise in the data are too poor to allow for any detailed profile fitting. Therefore, 
we generated a grid of photoionization models around the values suggested by our analysis 
of the LETG spectra (Netzer et al. 2002). Along with the models generated for the High 
Ionization component, the tables of ionic column densities predicted by Cloudy were used 
to construct a grid of high-resolution spectra, including the effects of bound-free absorption 
(Verner et al. 1996), and resonance and inner shell absorption lines (Behar & Netzer 2002). 
The spectral grids were used as an input to xspec in order to fit the XMM PN spectra. 


4.2. Spectral fits to XMM PN data 

Fig. 4 shows the 2001 April and November data relative to the powerlaw model (fit across 
the 3-5 keV band). This illustrates the strong curvature in both spectra. It is also clear that 
during November, the source was fainter and more heavily absorbed. We investigated the 
spectral shape for both epochs in light of insight from the Chandra spectral results. 

Prior to fitting we added to the model the set of emission features detected in our earlier 
analysis of the RGS data (Turner et al. 2003), which we suggested arose in gas outside the 
absorbers. Further, we discovered two emission features we had not picked up during our 
RGS analysis, one at ~ 0.83 keV, which is likely an unresolved blend of Fe L-shell and high 
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order Ovm resonance lines, and another at ~ 1.23 keV, which is most likely the Ne ix RRC, 
for which we measured a width of kT = 7.8 x 10 -2 keV. We used the xspec models zgauss 
and redge to parameterize these features in the model. We also included emission from 
the K-shell of Fe, the width was fixed to be a = 13 eV and energy fixed at 6.40 keV, as 
determined from HEG fits (see later). 

BeppoSAX data overlapped the XMM observations. The BeppoSAX data showed the 
spectrum of NGC 3516 to steepen to the canonical T ~ 2 in the 10-200 keV regime covered 
by the PDS spectra. This is indicative of the flatness being due to (unmodeled) complex ab- 
sorption or a high-energy cut-off to a flat powerlaw. In light of HETG detections of numerous 
absorption lines, we followed the former scenario. We found that when fitting BeppoSAX 
data simultaneously with PN data in xspec, the PN data almost entirely determined the fit 
parameters by statistical domination. Thus, to best estimate the strength of any reflection 
component relative to the primary spectrum, we used predominantly PN data from 3-10 
keV, tying the line strength to that of the reflector. In this way we obtained an estimate of 
i? ~ 1 for both epochs, representing reflection from steradians of material. This value 
was consistent with the spectral curvature observed in the BeppoSAX spectra. Therefore we 
included the reflected component fixed at R = 1 in all subsequent spectral fits. 

We constructed a model, using the parameters derived from the absorption lines for the 
High Ionization component and the Netzer et al. (2002) model for the UV absorber, compared 
to the April EPIC spectrum. The model drastically underpredicted the absorption at > 2 
keV and was unable to replicate the spectral curvature. Hence, there must be an additional 
large column of intervening gas (which we will refer to as the “Heavy” component). However, 
there did not seem to be any excess absorption in the soft band, which suggests that the 
Heavy component is in a high state of ionization (logU > 1) and that it may not fully cover 
the continuum source. 

Including this third (Heavy) component, we refit the April spectrum. In doing so, 
we assumed a priori that the UV component fully covers the X-ray continuum source, as 
suggested by Netzer et al. (2002) and Kraemer et al. (2002), as does the High Ionization 
component, as required by our absorption line fits. However, we let the covering factor of the 
Heavy absorber component be a free parameter. The radial velocity of the High Ionization 
component was fixed at -1100 km s -1 , as suggested by the observed energy of the H- and 
He-like absorption lines detected in the HETG data. The UV absorber was assumed to have 
an outflow velocity — 200 km s -1 (the mean from the strong UV components). We set the 
radial velocity of the Heavy component to —1100 km s -1 , however as we will show, the fit 
was not sensitive to the velocity. We obtained a statistically satisfactory fit with the complex 
model outlined above (Fig. 5a). The results are summarized in Table 2 (errors are x 2 -1-4.61). 



- 10 - 


Fig. 5b shows the effect of the Heavy absorber through a model plot with that component 
removed. 

Our spectral fitting returned a solution for Heavy that includes a high ionization pa- 
rameter and column density, but a covering factor of only 50%, which suggests this exists at 
a small radial distance, since the X-ray continuum source is < 1 lt-day in extent. The other 
consequence is that at such a low covering fraction, any absorption lines formed in Heavy 
will be undetectable since the continuum is dominated by the uncovered fraction. Note that 
the attenuation of the transmitted continuum by Heavy amplifies this effect. For example, 
at the energies of the high ionization lines 1.5-3 keV, the uncovered continuum contributes 
~ 75% of the total flux. Hence, we have no constraint on the velocity of this component, as 
no discreet absorption lines are detectable with HETG. 

We also tested the effects of screening of the continuum source by each of the compo- 
nents, as the photon input to each screen depends on any absorption of the nuclear flux 
by inner layers of gas. However, with such a low covering fraction we found little effect 
if the Heavy absorber was positioned closest to the nuclear source, wdth the UV and High 
Ionization components outside it. As this is a plausible radial ordering of the absorbers, we 
proceeded without worrying further about the effects of screening on the outer layers. 

The UV component of the model predicts column densities for H I, 1.4 xlO 17 cm -2 , 
C IV, 4.9xl0 16 cm -2 , and N v, 1.9 xlO 17 cm -2 , which are in agreement with the upper 
limits determined from the STIS spectrum of 2000 October 1, taken when NGC 3516 was in 
a slightly lower flux state (Kraemer et al. 2002). While Kraemer et al. derived a lower limit 
for S iv of 1.8 xlO 14 cm 2 , the UV component model predicts an ionic column density only 
~ 1/100 as large, which is consistent with the presence of additional lower ionization gas 
along the line-of-sight to the nucleus, which does not contribute appreciably to the X-ray 
opacity. The High Ionization component of the model predicts an H I column density of 1.4 
xlO 14 cm -2 , which may be detectable in Lya if the absorber covered the BLR, although the 
predicted C IV and N v were negligible. However, while there were components of Lya; at 
—1372 ±9 km s _1 and —994 ±16 km s -1 detected in the STIS spectrum (Kraemer et al. 
2002), the presence of C IV and N v lines at these velocities suggests that the High Ionization 
component is not seen in the UV, hence may not cover the BLR. The fact that the ratios 
of the ionic columns to that of H I were anomalously low in these components rules out the 
possibility that the Lya; lines arise in the X-ray absorber while the other lines form in lower 
ionization gas. 

The model of Netzer et al. (2002) explains the historic spectra of NGC 3516 utiliz- 
ing intrinsic absorption by the UV absorber, whose ionization changes in response to the 
continuum source. We modified the Netzer et al. (2002) model to account for the greater 
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complexity in absorption than was evident from previous datasets and, having achieved a 
good fit to 2001 April PN data, we fit the same model to the 2001 November PN data to 
determine the cause of the spectral variations. The November data were consistent with no 
change in column density for any absorption component, although there was a slight increase 
in covering factor (Table 2). All three absorption zones are consistent with a decrease in ion- 
ization parameter by a factor 1.5, which would be expected as the gas responds to changes in 
illuminating ionizing radiation. Fig. 6 shows the best-fitting models to April and November 
PN data, overlaid, illustrating the spectral change. 

Our final model, featuring three layers of intrinsic absorption, yielded a good fit to 
the overall shape of the X-ray spectrum. Another satisfactory aspect of the fit is that we 
recovered the canonical underlying photon index typical of AGN in the Seyfert class, and 
found a natural explanation for the F ~ 2 slope observed above 10 keV in PDS data. This 
gave us an additional level of confidence that we had correctly explained the flat observed 
spectrum of NGC 3516 in this low flux state. 


5. What of the Fe Ka Regime? 

5.1. The Line Core 

Having modeled the overall spectral shape we turned to a detailed examination of the 
residuals around 6.4 keV. From HEG data we found that comparison of the core of the Fe 
Ka line shows consistent fluxes in the 2001 April and November epochs, expected if the 
narrow line core is dominated by a contribution from the putative distant reflector (Fig. 7). 
HEG data yielded a core width a = lSl^eV, and energy 6.40 keV, as noted previously. The 
PN data, with superior signal-to-noise ratio, indicated an increase in core flux in November 
versus April, at > 95% confidence. If the core were simply velocity-broadened then a width 
of a ~ 13 eV yields a gas velocity 1400 km/sec. This would occur for gas lying ~ 100 lt-days 
from the black hole (assumed to have mass M ~ 2.3 x 10 5 * 7 M©, Wandel, Peterson, & Malkan 
1999), consistent with the suggestion of variability on timescales of months. However we 
note the mis-match of the XMM-based model (Fig. 7) with the HEG data, suggestive that 
the fit to the PN spectrum missed some of the core flux in 2001 April spectra. Thus the 
question of core variability remains open. 
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5.2. The Broad Diskline 

A long-standing question has been whether complex absorption of the continuum can 
mimic the broad red wing expected from Fe Ka contributions from the innermost accretion 
disk. It is clear that with our complex absorption model, whatever small excess of flux exists 
redward of 6.4 keV is not obviously characteristic of a diskline (Fig. 8a, b). 

We attempted to fit the April “excess” using the xspec model of a diskline following 
the prescription of Stella (1990) for a non-spinning black hole. This model produced an 
improvement d% 2 = 31 with inner and outer radii fixed at 6 and 1000 gravitational radii, 
respectively, disk inclination fixed at 30° and emissivity index at q = —2.5. The line energy 
was 6.40 < 6.45 keV, equivalent width 4311*72 eV. Thus we conclude the data are consistent 
with the presence of a substantial diskline, although the spectra do not provide compelling 
evidence for such. 


5.3. Red or blue-shifted narrow components of Fe Ka emission 

Turning to the question of narrow line emission, we examined the spectra for narrow 
red or blueshifted Fe ka components, as reported previously for this source (Turner et al. 
2002, Bianchi et al. 2004). We split the PN data from 2001 April, dividing the spectra ~ 30 
ksec into the 49 ksec exposure, splitting the good interval into two sections, between which 
the source flux increased by ~ 12%. 

Fig. 8a and b show the Fe Ka profile from the two time-selections sampling changes in 
flux within each epoch of data. There is evidence for a very weak feature at rest-energy 6.08 
keV, which was reported by Bianchi et al. (2004). As previously reported for November data 
(Turner et al. 2002), spectral variability is evident in the Fe Ka regime down to timescales 
of a few tens of ksec and these are possibly related to changes in continuum flux. While it 
is tempting to attribute the variable spectral shape in the Fe Ka regime (most evident in 
November) to changes in a broadened line from the accretion disk, our insight from HEG 
spectra leads us to attribute this variability to rapid (tens of ksec or less) changes in the 
energy-shifted narrow components of Fe Ka. We also note the presence of additional varying 
features not previously reported, e.g. close to 3 keV (Fig. 8b). While it is tempting to 
speculate on their origin, much longer XMM observations are needed to progress on reliably 
characterizing the line identifications and behaviors. 
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6. Discussion and Conclusions 
6.1. Nature of the Absorbers 

I 

We find a model for NGC 3516 which has three intrinsic layers of absorption with sig- 
nificantly different columns and ionization-states, and, when determined, outflow velocities. 
The UV absorber has previously been detected in both UV and soft X-ray regimes. A High 
i Ionization absorber is required to produce the strong H- and He-like Mg, Si, and S lines 

detected in the 2001 April HETG spectrum. Finally, the Heavy absorber has been detected 
previously but we find the first constraint on its covering fraction. 

| 

The first evidence for the Heavy component came with the detection of an iron edge 
at ~ 7.8 keV in Ginga data (Nandra k Pounds 1994) taken when NGC 3516 was observed 
several times in flux-states close to that found for 2001 April. Those authors found a column 
of a few x 10 23 cm -2 for that gas, with similarly high ionization-state to the Heavy component 
in our models. Costantini et al (2000) also found the presence of a large column of highly 
ionized gas in BeppoSAX spectra, although at Ng ~ 2 x 10 22 cm~ 2 this is different to the 
i Nandra & Pounds (1994) solution and resembles our High Ionization rather than our Heavy 

component. Since the High Ionization component would be undetectable at low resolution 
if a significantly larger column of gas were present, the discrepancy may indicate that the 
Heavy column was absent during the earlier BeppoSax observations. 

Based on our models, we can determine the densities, tih, of the High Ionization and 
Heavy components as a function of their radial distance, r. For the High Ionization absorber, 
log£7 = 1.7, hence for a luminosity in ionizing photons of ~ 2 x 10 53 s -1 (Kraemer et 
al. 2002), rn 1 ^ 2 ~ 1.0 x lO^cm -1 / 2 , while for the Heavy absorber, log U =1.2 and rn# 2 ~ 
1.8 x 10 20 cm -1 / 2 . Although we cannot deconvolve the densities and distances of the absorbers, 
it is implausible that the physical depth of the absorbers (Ar = n# /N h) could be more than 
a small fraction of their radial distance. For example, for the High Ionization component, 
if A r/r = 0.1 then n# = 2.5 x 10 6 cm -3 and r = 6.3 x 10 16 cm, which is approximately 
equal to the lower limit derived for the UV absorber (Kraemer et al. 2002). The same line 
of reasoning suggests that the Heavy component lies somewhat closer to the central black 
hole, which is consistent with partial covering. In fact, since the size of the X-ray continuum 
source is < light-day, one might expect the radial distance of the Heavy absorber to be of 
the same order of magnitude. Alternatively, partial covering may result if the absorber is 
patchy or filamentary. 

Although we cannot constrain the densities of the Heavy and UV components, the ratios 
of their mean electron temperatures to their ionization parameters can be used to estimate 
their relative thermal pressure. The model predicts T = 10 5,63 K, 10 5 " K and 10 4,44 K for 
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the Heavy, High Ionization, and UV absorbers, respectively. The thermal pressure of the UV 
absorber is ~ 6 times that of the Heavy absorber and 9 times that of the High Ionization 
absorber. Hence, these three components cannot be in thermal equilibrium with respect 
to each other. To examine the thermal stability of the absorbers, in Fig. 9 we show the 
relationship of T and the so-called “pressure” ionization parameter U /T (generally known 
as an S-Curve). As discussed most recently in Krolik & Kriss (2001), there are two regions 
where a photoionized gas is stable to thermal perturbations: at low ionization/temperature, 
when line cooling is efficient, and high ionization/temperature, when thermal balance is 
achieved via Compton processes. In the intermediate region, the gas can be unstable if 
dT/d(U/T) < 0. There is no such region for our models, in part due to the additional 
cooling provided by the over-abundance of nitrogen (in contrast, see the S-curve in Kraemer 
et al. 2002). Nevertheless, as a result of the large changes in temperature which result from 
small changes in U, the intermediate region is only marginally thermally stable (for a more 
detailed discussion, see Krolik & Kriss 2001). While the UV absorber is within the line- 
cooled region, the two highly ionized components lie in the intermediate region, and should 
undergo large changes in temperature in response to variations in the incident continuum. 

Given the high ionization state of the High Ionization and Heavy absorbers, they are 
quite transparent to the incident continuum, with the main source of opacity being bound- 
free transitions from highly ionized O, Ne, Mg, Si, S, and Fe. Our models predict an average 
force multiplier, the ratio of the radiation pressure force from all absorption and scattering 
process to that caused by Compton scattering, ~ 3. As a consequence, the central source 
would have to be radiating at roughly 30% of its Eddington Luminosity in order to generate 
a radiatively driven outflow. Also, our models predict mean temperatures which are too 
low for a thermal wind to achieve such a high velocity (Balsara & Krolik 1993). This 
suggests the possibility that these absorbers are associated with a hydromagnetic (MHD) 
wind (e.g., Emmering, Blandford, &; Shlosman 1992). As noted above, their high ionization 
state strongly suggests an origin close to the continuum source and Emmering, Blandford, 
& Shlosman (1992) predict that the outflow velocity in an MHD flow should increase as the 
radial distance of the launch pad decreases. Furthermore, based on their models of MHD 
flows in the BLR gas in NGC 5548, Bottorff et al. (1997a) predict line-of-sight velocities 
of the same order as that observed for the High Ionization absorber. Finally, as discussed 
above, these absorbers are only marginally thermally stable, although they appear to have 
been present in earlier observations. If the absorbers were permeated by magnetic fields of 
sufficient strength (i.e., on the same order as their thermal pressure), the magnetic pressure 
would act to stabilize them (Bottorff, Korista, & Shlosman 2000). The presence of internal 
magnetic fields would also suggest an MHD origin. 

In principal, given the line-of-sight distance of the absorber from the black hole, the 
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line-of-sight velocity of the material, the disk orientation with respect to the observer, and 
the black hole mass, it may be possible to constrain parameters of a Blandford & Payne 
(1982) type MHD disk wind model applied to NGC 3516. This is the subject of future 
research efforts. We note that an early attempt to constrain Blandford & Payne (1982) 
model parameters in NGC 5548, using its BLR Civ emitting region, is described in Bottorff 
et al. (1997b) and was relatively successful in describing emission line shape and variability. 
This model was extended to the UV/X-Ray warm absorber of NGC 5548 but analysis was 
hampered by the then relatively low resolution X-Ray data of the pr e- Chandra, pr e-XMM 
era (Bottorff, Korista, & Shlosman 2000). We also note that the masses of the black holes in 
NGC 3516 and NGC 5548 are similar. It would therefore be tempting to directly apply the 
results of the NGC 5548 model to NGC 3516. We caution however that slight differences 
in Blandford & Payne (1982) MHD flow parameters and/or orientation can lead to large 
differences in the radial velocity and line-of-sight position of absorbing material. 


6.2. The Fe Emission 

We note with interest the remarkable similary between the Fe Ka profile of NGC 3516 
and that of Mrk 766 (Pounds et al. 2003a, Turner et al. 2003). Mrk 766 and NGC 3516 are 
two of very few Seyfert galaxies known to exhibit narrow and highly-redshifted Fe emission 
lines (Turner et al. 2002, Turner et al. 2003). Another similarity is that Mrk 776 also 
shows evidence for a large column of ionized gas (Pounds et al. 2003a) with N H ~ 2 x 
10 23 atom cm -2 . If the emission is from ionized Fe then it would be natural to expect the 
heavy absorption and emission from ionized Fe to be observationally linked. It may be that 
we can only see the weak narrow Fe lines in sources with the largest columns of circumnuclear 
gas. If the narrow redshifted Fe emission lines (Turner et al. 2002) do originate in the highly- 
ionized gas components then previous speculations will have underestimated their velocity. 
The emission from the K-shell of Fe in NGC 3516 is expected to peak at Fe xx (~ 6.51 
keV), perhaps raising estimates of the outflow velocity for the emitting gas responsible for 
the shifted features reported by Turner et al. (2002) up to ~ 0.16c, much higher than outflow 
velocities estimated for the absorption lines reported here. 

Other AGN have shown evidence for zones of high-ionization gas, e.g. NGC 3783 
(Reeves et al. 2004) which shows an Fe resonant absorption line, where again, the inclusion 
of the absorber significantly reduces the inferred equivalent width of any broad component 
to the diskline. Even higher column (~ 10 24 atom cm -2 ) outflows of gas have also been 
suggested with outflow velocities up to 0.1c in QSOs such as PG 1211+143 (Pounds et al. 
2003b), PDS 456 (Reeves et al. 2003) and the NLSyl NGC 4051 (Pounds et al. 2004). In 
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those cases the column density is high enough that absorption edges + Fe K-shell resonant 
absorption lines can be detected. 

While current data leave us tantalizingly close to answering the question as to whether 
high-column, high-ionization absorbers are masquerading as broad disklines, a conclusive 
answer will have to wait until long exposures are taken on AGN using the microcalorimeter 
on ASTROE-2. 
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Table 1. X-ray Flux 


Epoch 

2-10 keV a 

0.5-10 keV a 

7 keV b 

Observed fluxes 

2000 Oct 

1.01 

0.17 

0.78 

2001 April 

2.22 

0.43 

1.72 

2001 Nov 

1.55 

2.70 

1.35 

Absorption-Corrected Fluxes 

2000 Oct 

1.27 

0.70 

0.82 

2001 April 

2.90 

1.81 

1.82 

2001 Nov 

2.11 

1.11 

1.45 

Absorption-Corrected PowerLaw Flux 

2000 Oct 

1.04 


0.61 

2001 April 

2.71 

1.78 

1.62 

2001 Nov 

1.80 


1.15 


a Fluxes as ergs cm 2 s 1 in units of 10 11 . 

b Fluxes as keV cm -2 s ~ l keV~ l in units of 
IQ- 3 . 
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Table 2. PN fit parameters 


2001 April 2001 Nov 


r 

Cov a 
Log N h,uv 
Log U xjv 
Log Nh,m 
Log U hi 

Log N h , heavy 
Log U heavy 

E FeK (keV) 

nFeK C 

X 2 /d.o.f. 


1.82±0.01 

44±6% 

21.85±0.15 

-0.80i0.02 

22.2/ b 

1-7/ 

23.40±Hs 

i oi +0.02 

i * Zi -0.04 

6.40/ 

3.64l°i« 

1801/1521 


1.77i0.02 
58±5% 
21.60i0.15 
— 0.99±0.04 
22 . 2 / 
1.5/ 

23.40±0.04 

1 . 091^4 

6.40/ 

A Q1 +0.22 
— 0.36 

3248/2770 


a Fraction of nuclear source covered by the 
heavy absorber 

b / denotes that a parameter was fixed in 
the fit 

c Photons in the line in units 10 -5 
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0.65 0.7 

Rest Wavelength (nm) 



Fig. 1. — The panels show key sections of HETG data, where we have combined HEG and 
MEG data from the ±l st -order spectra. The most prominent absorption lines are evident 
in the different panels, with a) Sxvi; b) Sixiv and c) Mgxn. The dotted lines show where 
those and other absorption lines would be expected. We have ignored some artifacts of the 
instrument which mimic absorption features; for example, the Mg IX region falls on the 
S3/S4 chip gap in the positive first order spectrum, producing a depression in the observed 
spectrum. 




Observed Redshift 
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Fig. 2. — Velocity dispersion of the absorber versus outflow velocity. The 90% confidence 
contours are shown for the kinematic parameters of the lines from the HETG data. The 
absorption lines are blueshifted relative to the systemic redshift, with Sxvi showing a higher 
blueshift than the other lines, as discussed in the text. It is also evident that the lines yield 
no significant constraint on the velocity dispersion of the absorbing gas, and so a value of 
300 km s -1 has been adopted. The FWHM for the MEG and HEG are shown for reference 
for the lowest and highest wavelengths (these are for Sxvi and MgXli lines). 


-21 - 



-10 0 10 


Offset 

from Systemic 

(0.01 A bins) 


Fig. 3. — Velocity spectra, binned at O.OlA for the three strongest lines detected in the 
2001 April MEG (only) spectra from obsid 2080. The red line shows the model fit to the 
absorption feature convolved with the MEG spectral resolution and assuming a velocity 
dispersion of 300 km s -1 . 
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channel energy (keV) 


Fig. 4. — Ratio of PN data from 2001 April (top line, black) and November (bottom line, 
red) to a model which represents the mean power-law continuum of the two epochs fit over 
in the 3-10 keV regime. This plot illustrates the flattening of the spectrum in November 
compared to April, due to the drop in ionization-state of the absorber as the source flux 
decreases. 
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Fig. 6. — The unfolded spectrum (from fits to the EPIC PN data) compared to the model 
discussed in sections 4 and 5. The green line (top) represents the 2001 April data while the 
red line (bottom) is based upon 2001 November data. 







Observed Energy (keV) 


Observed Energy (keV) 


Fig. 8. — Ratio of data/model for our best-fitting model from Table 2. Each of the PN 
datasets have been split in time to sample different flux states, as discussed in the text. The 
red line represents the higher flux state in each case. Rapid changes evident in the profile 
during November are due to rapid variations in the flux of narrow and redshifted Fe lines 
(Turner et al. 2002). April data do not allow us to sample such a wide range of flux and 
while we see evidence for a previously-reported line at ~ 6.1 keV (Bianchi et al. 2004) no 
significant variability is evident within the good section of PN data. 



Fig. 9. — Thermal stability curve for models generated with our assumed SED. The flat 
sections, at low and high temperature, are the stable line-cooled and Compton-cooled regions, 
respectively. Note that, while there are no thermally unstable (dT / d(U /T) < 0) regions, in 
the region where the curve is nearly vertical small changes in ionization can cause large 
changes in temperature. The locations of our three models for the 2001 April spectra are 
shown. 
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